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ABSTRACT 

We report on the first detection of a quasi-simultaneous radio-X-ray flare of Cygnus X-l. The detection was 
made on 2005 April 16 with pointed observations by the Rossi X-ray Timing Explorer and the Ryle telescope, 
during a phase where the black hole candidate was close to a transition from the its soft into its hard state. The 
radio flare lagged the X-rays by ~7 minutes, peaking at 3:20 hours barycentric time (TDB 2453476.63864). 
We discuss this lag in the context of models explaining such flaring events as the ejection of electron bubbles 
emitting synchrotron radiation. 

Subject headings: X-rays: stars - X-rays: binaries - Black Hole Physics - accretion, accretion disks 



1. INTRODUCTION 

With the increased availability of simultaneous radio and 
observations in the last decade, there is now a large amount of 
evidence available pointing towards a very close interaction 
between the accretion disk and the jet in black hole X-ray bi- 
naries and active galactic nuclei (AGN). Most convincingly, 
this disk-jet interaction has been shown for microquasars, 
i.e., black hole binaries with strongly relativistic jets such as 
GRS 1915+105 or GRO J 1655 -40. In these systems, the cor- 
related flaring in the X-rays, optical/infra-red, and radio seen 
at certain times is generally interpreted as the evidence for 
(ballistic) ejection events of synchrotron radiation-emittin 
electron bubbles (iRothstein et al.1 120051: iFender & B ellon: 
12004 IKlein-Wolt et al.1 l2002t lEikenberrv et al.1 1 1 9981 and 
therein). In this model, the X-ray flare represents the ejection 
of the synchrotron radiation emitting bubble, which then adia- 
batically expands within the jet flow and cools down, resulting 
in the pea k of the emission shifting downwards in frequenc y 
with time (van der LaanlfT%6l iHjellming & Johnston| [l988). 
Simultaneous broadband observations of such events, which 
show minute-long delays betwe en the different wave bands , 
are consistent with this picture (iMirabel & Rodriguez! 1 19941: 
Mirab el et alJI 19981: iPoolev & Fendeill 19971 lEikenberr v et al.l 



1998b . The model has also been confirmed by proper motion 



measurements in the radio, which reveal intrinsic jet sp eeds of 
>0.57cforGRS 1915+105 (iMiller- Jones et al J 120051) . Com- 
parable behavior was also detected in 3C 120, suggesting that 
similar ejections also occur in active galactic nuclei, on corre- 
spondingly longer timescales (Marscher et al. 20Q2|). 
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For black hole binaries with weakly relativistic jets, the ev- 
idence for jet-disk-interaction is less direct. This evidence 
includes the correlation between X-ray states an d radio emis- 
sion i n black hole transien ts (e.g., in GX 339-4 Corb el et"aLl 
l200l iBelloni et all l2005h and the success of modeling 
the radio to X-ray broad band spectrum of black hole 
candi d ates with outflow-domi nated models dMarkoff et all 
120051: iMarkoff&Nowakl 12004 and therein). Further- 
more, at least for Cygnus X-l, there is also evidence for 
the presence of an e nergetically s i gnific a nt, strong outflow 
(Stirli ng et al.1 \l99§. IStirling et al] I200TI iGallo et alj 12001 
IMiller- Jones et alj|2006l) . A relativistic jet with v > 0.3c has 
been associated with radio flares in this system (Fend eret alj 
l2006l) . 

Apart from GRS 1915+105, however, none of these obser- 
vations shows direct evidence for a causal connection between 
the X-rays and jet on time scales of minutes. Prompted by 
this lack of quasi-simultaneous short-term radio-X-ray corre- 
lations, in 1998 we initiated a long term monitoring campaign 
of Cyg X-l with the Rossi X-ray Timing Explorer (RXTE) 
and the Ryle telescope. Biweekly 3-10ksec long simultane- 
ous observations started in 1999. Previous searches for flares 
in campaign data taken between 1999 and mid-2003 did not 
reveal evidence for coherent short term activity in both bands, 
although a significant correlation on t ime scales of weeks 
was found, especi ally above ~10keV (Glei ssner et alj[2004t 
Wilms 'etaT]l2006l) . In this Letter, we report on the observa- 
tion made on 2005 April 16, in which the first clear quasi- 
simultaneous radio-X-ray flare was detected in Cyg X-l. The 
remainder of this Letter is structured as follows. In §|2]we de- 
scribe the observations, followed by the analysis of the flare 
in §[3] We discuss the physics of the flare in the context of 
emission models for the radio and X-ray emission in $4] 

2. OBSERVATIONS AND DATA REDUCTION 

We use data from both instruments on-board t he RXTE, the 
low e nergy Proportional Counter Array (PCA; Jaho da et al.1 
120051) an d the High Energy X -ray Timing Experiment 
(HEXTE; iRothschild eTail Il998l) . The data analysis was 
performed using the standard RXTE data analysis software, 
HEASOFT 6. 1 .2. Spect ral fitting was performed with XSPEC 
11.3.2aa dArnaudlll996l) . 

A crucial part of the observation happened during the early 
phase of the RXTE observation, shortly after the source rose 
above the Earth's horizon. Due to auroral emission in the far 



2 



Wilms et al. 



25(1 



200 



JD-2450000 

3200 3300 3400 



3500 



0.63 



JD-2453476 

0.64 0.65 0.66 



0.67 




JASONDJFMAMJ 

2004 2005 

Fig. 1 . — RXTE-ASM (lower curve, left y-axis) and Ryle telescope 1 5 GHz 
(upper curve, right y-axis) lightcurves of Cyg X-l from 2004 March until 
2005 June. The time of the simultaneous observation is indicated by the 
dashed vertical line. Gaps are shown if their duration is >4d. 

ultra-violet and soft X-rays and due to cosmic ray reprocess- 
ing in the hard X-rays, the Earth's atmosphere is not com- 
pletely X-ray dark. The typical 2.5-20 keV X-ray flux at 
the typical magnetic latitude of the RX TE orbit is too low , 
however, to influenc e our measurements dPetrinec et alJ2000t 
ISazonov et al.1120 07). We therefore use all data taken whilst 
the source was > 1 above the Earth's horizon and had a source 
offset of <(f.0l. We use PCA data from the standardl mode, 
which gives the full 2.5-128 keV PCA count rate with a time 
resolution of 0.125 s and no energy information, and from 
the standard2f mode, a binned data mode with a 128 chan- 
nel energy resolution and a time resolution of 16 s. X-ray 
light curves were extracted with the intrinsic time resolution 
of each mode and then barycentered and rebinned. 

The Ryle telescope data were taken at 15 GHz with a 
time resolution of ~8s. The typical lcr uncertainty of the 
radio measurements is 9mJy. The observations are inter- 
rupted every ~ 1600 s for phase calibration observations of 
J2007+4029. The amplitude c alibration o f the R yle data cor- 
responds to the flux scale of Baa rs et ail (1 19771) and is per- 
formed using nearby ob servations of 3C48 and 3C286. See 
iPoolev & Fender! (1 19971) for further information on the Ryle 
telescope. 

3. A QUASI-SIMULTANEOUS RADIO-X-RAY FLARE 

As shown in Fig. Q] 2005 April marks the possible end 
of a longer X-ray flaring episo de of Cyg X- 1 that started in 
early 2004 (I Wilms et all 120061) . While clearly defined ra- 
dio fl ares are not uncommon in Cyg X-l (e.g., iHjellmingl 
1 1973b . increased radio emission and radio flaring are gener- 
ally seen when the source is in the intermediate state between 
the hard and the soft states, while the r adio is weak once 
the X-ray source approaches the soft state ( Wilms Tt al.ll2.Q06l 
and therein). At the time of our pointed observations, the soft 
X-ray flux had just come down from a large flare. Shortly 
after the observation, the Id averaged 15 GHz flux peaked, 
reaching a maximum of ~30 mJy, close to the brightest radio 
flux of Cyg X- 1 during 2004/2005 . 

Figure shows the 15 GHz radio flux and the RXTE PCA 
count rate lightcurve measured on 2005 April 16. Close to the 
start of the observation, a radio flare is readily apparent. The 
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Fig. 2. — Barycentric 15 GHz flux density (left axis, dotted line) and 2— 
128 keV count rate lightcurves (right axis, solid line) of Cyg X-l for the flare 
event of 2005 April 16, both binned to a resolution of 32 s. The bottom axis 
gives the barycentric time in hours, with tick-marks separated by 10 minutes. 

total duration of the flare is ~15 minutes. During this interval 
the 15 GHz flux increased by a factor of ~3 to a peak radio 
flux of 70 mJy. This radio flux is among the highest seen dur- 
ing the Ryle monitoring 8 . Previous radio flar es, however, did 
not occur during pointed RXTE observations (Gleissn er et al.1 
120041: IPender et al.l2006l) . and the source monitoring provided 
by the RXTE-ASM is too coarse to pick up such short lived 
X-ray events. 

RXTE started observing Cyg X-l about 10 minutes before 
the peak radio flux. The X-ray lightcurve shown in Fig. [2] 
shows a similar shape to the radio one, although with more 
substructure. The earlier maximum of the X-ray flare did 
not allow RXTE to catch the start of the X-ray flare, or de- 
termine whether the maximum X-ray flux seen is indeed the 
peak of the X-ray flar e. A cro s s corr elation (CCF) analysis 
using the algorithm of IScargla d 1989b reveals a413 + 165s 
time lag of the radio with respect to the X-rays, where the 
lcr uncertainty was determined using a standard bootstrap- 
ping method with 1000 realizations. Other approa ches to cal- 
culate the CCF for non-unifor mly sampled data d Alexander! 
1997; Edel son & Kro lik 1988) give essentially the same re- 
sult. With a maximum lScarglel (ll989h -CCF of 0.38, this anal- 
ysis formally confirms the general similarity of the X-ray and 
radio lightcurves (Fig. [3]). Since the substructure of the X-ray 
lightcurve, i.e., the two smaller flares after the main flare, is 
clearly different from that in the radio flare and since the start 
of the X-ray flare is not covered by our observations, the peak 
CCF value is not higher. For the same reasons, the formal 
uncertainty of the lag measurement is rather large. 

To characterize the shape of the radio flare, we fit the radio 
data (rebinned to a resolution of 8 s) with the sum of a linear 
flux trend and a Gaussian representing the flare, 



f{t) - a(t-t') + b+A exp - 



(t-TY 



(1) 



where t' is a reference time, taken as the center of the time 
interval analyzed (f = JD 2453476.64635). The barycen- 
tric time of maximum flux occured at T — JD 2453476.63864 
with an uncertainty of 15 s (uncertainties are at the lcr level 
for one interesting parameter). The peak flux of the flare 

8 The most exceptional radio flare was that of 2004 February 20, which 
reached a peak flux of 140 mJy at 15 GHz, the larg est flux ever seen for this 
source with the Ryle telescope I Fender et al. 2006). 
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Fig. 3. — Cross-correlation function (Scargle 1989) for the radio data with 
respect to the RXTE -PCA data. A negative time lag indicates the radio lag- 
ging the X-rays. 

component is A = 39.2 + 2.5 mJy and the width of the com- 
ponent is <x = 210 + 16 s. The flare is superimposed to a 
continuum with b — 25.3 + 0.7 mJy, decreasing linearly with 
a = -13+2 mJy h 1 . The high quality of the fit is indicated by 
its low reduced ;f 2 (y 2 . = 1.11 for 233 degrees of freedom). 
The radio flare is therefore symmetric around its maximum. 
This symmetry is a marked differ ence compared to th e asym- 
metric flare of 2004 February 20 dFender et al.ll2006l) . 

Modeling the RXTE-PCA light curve is complicated by 
the flare already being in progress when the measurements 
started. Furthermore, contrary to the radio data, where the 
scatter in the lightcurve is mainly due to the measurement 
uncertainty, the X-ray data are dominated by strong low- 
frequency noise onto which the X-ray flare is superimposed. 
Consequently, the empirical model of Eq. (Q~|) does not result 
in a good description of the X-ray data. 

To study the spectral evolution of Cyg X-l during the 
flare, we perform a spectral analysis of the 2.5-20 keV 
PCA standard2f data at 16 s time resolution using a sim- 
ple photoabsorbed powerlaw, which proves sufficient to de- 
scribe the spectrum at this lower signal to noise level. For 
spectra taken during the flare, at PCA count rates above 
1000 counts s _1 PCLT 1 , the mean power law index F = 2.10 ± 
0.03. The spectrum hardens outside of the flare to T = 
1.98 ± 0.03 (errors given are the standard deviation of the 
individual power law fits to the standard2f spect ra), a value 
typica l for the intermediate state of this source dWilms et alJ 
120061) . 

That Cyg X-l was in the intermediate state on the day 
of the flare can also be confirmed by modeling the 2.5- 
150keV PCA and HEXTE spectrum of an RXTE observa- 
tion performed 4h after the flare (to avoid possible "con- 
tamination" by the flaring activity) with the sum of a pho- 
toabsorbed (A^h = 6 x 10 21 cm 2 , held fixed), exponentially 
cut-off broken power law model. This empirical model has 
been shown to give a good charac terization of the spectral 
shape of Cyg X-l (Wilms et al. 2006). The spectral parame- 
ters are a lower photon index Ti = 2.01 ± 0.01, breaking at 
£break = 10.0+°^ keV into a power law with T 2 = At 
E cut = 26 + 3 keV the exponential cutoff starts with a folding 
energy of EfM = 1 37 ± 1 1 ke V (all uncertainties are at the 90% 
level). In addition, a Fe Ka line from neutral iron is present 
with an equivalent width of 135 eV. The 3-10keV source flux 
is 6.17 x 10~ 9 ergscrrT 2 s _1 . The parameters of the contin- 
uum are again consistent with an intermediate state and fits 
well with the empirical picture that radio flaring in black hole 
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Fig. 4. — 1-128 keV RXTE-PCA count rate versus hardness ratio for Cyg 
X-l observations from the 2000-2005 RXTE monitoring (black) and for data 
taken during the 2005 April 16 flare (red). The dotted line indicates the lo- 
cation at which observations with maximum radio flux are found in the long 
term monitoring. The hardness ratio was determined from background sub- 
tracted PCA top layer data from channels 8-14 and 15-25, respectively. A 
color version of this figure is available in the online edition of ApJL. 

candidates occurs most frequently in this state ( Fender et al.1 
2004: lWilmseTaIll20 06). 

To allow the interpretation of the observed softening during 
the flare with the general behavior of Cyg X-l, Fig. [4] shows 
the X-ray hardness intensity diagram for the PCA top layer 
standard2f 16 s spectra in the context of the pointed RXTE 
observations of the monitoring campaign. For black hole tran- 
sients this diagram is seen to have an approximate q-shape 
dBelloni et a l. 2006, and therein). As a persistent hard state 
source, Cyg X-l is typically found in the top right corner of 
the diagram. Outside of the flare, the source is situated at a 
hardness of ~0.55 with a typical 1-128 keV PCA count rate 
of ~900 counts s 1 PCLT 1 . During the flare the source soft- 
ens and brightens. It leaves the region of the diagram where 
Cyg X-l is usually found during our monitoring campaign, by 
moving to higher count rates, for this hardness, than usually 
observed. 

Flaring behavior in Cyg X-l is usually observed when- 
ever the source is close to the "jet line" in its hardness in- 
tensity diagram, while the radio flux gets quenched once the 
source moves away from the jet line to th e left of the dia- 
gram (tGalioetiDllQOl lend^^^^S. In spectral fits 
based on the eqpair-model of lCoppil £1999), the radio fluxes 
of Cyg X-l are at their maximum when th e compactness ra- 
tio 4/4 ~ 3 (Fig. 16 of lWilms etaill2006l ; the compactness 
ratio is a measure for the relative importance of the energy 
dumped into the Comptonizing plasma and that dissipated in 
the accretion disk), corresponding to a soft power law index of 
Ti ~ 2.1. From our database of spectra of Cyg X-l, we find 
that these observations have a (7.1-1 1.7 keV)/(4. 6-7.1 keV) 
hardness ratio of 0.52, indicated by the dotted line in Fig. [4] 
During the 2005 April 16 observation, the source was there- 
fore close to this line of maximum radio flux and approached 
it asymptotically during the flare. 

4. DISCUSSION 

In this Letter we have presented the first evidence for a di- 
rect relationship of the X-ray and radio emission in Cyg X-l 
on timescales of minutes. The data show the radio to lag the 
X-rays by 413 + 165 s and the X-ray spectral shape to ap- 
proach the X-ray hardness ratio where the source is typically 
found at its largest radio flux in our long term monitoring, 
hinting towards a general similarity of the physics of individ- 
ual flare events and the overall radio-X-ray connection. Al- 
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though the X-ray data do not cover the start of the X-ray flare, 
explaining the rather large uncertainty of the lag determina- 
tion, the morphological similarities between the X-ray and ra- 
dio lightcurves also suggest that the same event is observed in 
both wavebands. 

Similar events in microquasars show that lags with 
timescales of several 100 s are typical for the coherent be- 
havior of these systems, s uch as a lag of 310+20 s between 
the soft X-rays and the IR dEikenberrv et alJI 19981) and 80 s 
between the X-rays and the radio dPoolev & Fenderll 19971) in 
GRS 1915+105. The timescale observed in Cyg X-l allows 
us to place an upper limit to the physical separation of the X- 
ray and radio emitting regions of the accretion/ejection flow. 
We assume that the emission coincides with the imaged jet 
and that the jet is perpendicular to the orbital plane of the 
HDE 2 26868/Cyg X-l sy stem (although this is not a priori 
certain; | Maccarong |2002|). which has an inclination of 30° 
dGies & Boltonlll982hlDolanlll992l) . Taking light travel time 
effects into account, for jet speeds of 0.3 c, the lower limit im- 
plied by observation s of the transient jet ejection discussed by 
iFender et al.l (120061) . the measured delay implies a separation 
of 1 . 1 + 0.5 AU between the location of X-ray and radio emis- 
sion. If the jet is relativistic instead, with a speed of 0.99c, 
the distance increases to 5.8 + 2.3 AU. Note that similar val- 
ues for the length of the jet are obtained by considering that 
the ~ 10 minute duration of the radio flare is roughly equal 
to the dynamical timescale of the jet. Assuming a distance 
of 2.5 kpc, these values imply to a maximum projected angu- 
lar separation between the X-ray and radio emitting region of 
10~ 3 mas. 

What is the phy sics of the observed event? In the model of 
IFender et al.l (|2004) for transient radio events, the inner edge 
of a thin accretion disk is posited to move rapidly towards 
the black hole. The temperature at the inner edge therefore 
increases, leading to a softening of the source in the RXTE- 
PCA as more disk photons enter the instrument's band pass. 



This X-ray flare is then followed by the ejection of an elec- 
tron bubble, which rapidly expands, producing the observed 
radio emission. At least qualitatively, this behavior and also 
the time scales deduced above seem to agree with our obser- 
vations, although the model was originally invented for the 
large scale variability of black hole candidates and not for 
such short events as the one discussed here. Note that while 
the main flare dominates the measured time lag and therefore 
the sizes estimated above, it is followed by two short spikes, 
which are both present in the radio and the X-ray lightcurves, 
but at different time delays. These spikes could indicate that 
more than one blob of material was ejected at different speeds 
but cannot be separated once the blobs have expanded and 
their radiation peaks in the radio. Such a behavior could be 
typical for flares in Cyg X-l, since the large radio flare of 
Cyg X- 1 from 2005 Feb 20 also shows very little substructure 
dFender et al.ll2006l) . 

The lack of further detections of radio-X-ray flares in over 
1.5 Msec of simultaneous radio-X-ray data precludes a more 
detailed discussion of the properties of the ejected material 
from comparing the source behavior during different flares. 
The observation of the flare itself, however, stresses the im- 
portance of long-term multi-wavelength campaigns to detect 
such rare events which are necessary to further our insight into 
the physics of the emission from black holes. 
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